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INTRODUCTION 


High  temperature  polymer  composites  are  receiving  special  attention 
because  of  their  potential  applications  to  high  speed  transport  airframe  structures 
and  aircraft  engine  components  exposed  to  elevated  temperatures.  Typical 
composites  considered  for  these  applications  are  carbon/polyimide  composites 
such  as  IM7/K3B.  The  long  term  durability  of  these  materials  is  of  great  concern. 
It  is  related  to  and  affected  by  physical  aging  of  the  polymeric  matrix,  viscoelastic 
deformation  (creep)  and  hygrothermomechanical  degradation.  It  is  therefore 
important  to  investigate  the  processes  involved  and  any  accelerated  test 
methodologies  that  will  yield  valid  predictions  of  long  term  behavior. 

The  objective  of  this  investigation  is  to  characterize  the  constitutive 
behavior,  damage  mechanisms  and  damage  development  in  high  temperature 
polymer  matrix  composites  and  develop  models  for  long  term  behavior  under 
various  environmental  and  loading  conditions.  The  approach  followed  is  as 
follows: 

1 .  Select  candidate  polymer  composite  materials. 

2.  Characterize  unidirectional  lamina. 

3.  Perform  aging  and  viscoelastic  characterization  under  isothermal  and 
nonisothermal  conditions. 

4.  Study  static  and  fatigue  behavior  of  laminae  and  laminates  at  various 
cyclic  stress  amplitudes  and  temperature  levels. 

5.  Study  failure  modes  and  damage  development. 

6.  Determine  stiffness  degradation  as  a  function  of  stress  amplitude,  fatigue 
cycles  and  temperature. 

7.  Develop  accelerated  testing  methodology  and  life  prediction  models. 

Three  composite  material  systems  were  selected  for  the  investigation: 
IM7/K3B  carbon/polyimide,  IM7/RP46  carbon/polyimide  and  IM7/977-3 
carbon/epoxy.  These  systems  were  characterized  under  monotonic  loading 


conditions  at  room  temperature.  Creep  and  physical  aging  of  one  of  these  systems 
was  investigated.  Processing,  characterization  and  residual  stresses  as  well  as 
failure  mechanisms  and  damage  development  were  also  investigated  under 
monotonic  and  cyclic  loading  conditions. 

BASIC  MATERIAL  CHARACTERIZATION 
Carbon/Polyimide  (IM7/K3B) 

This  material  consists  of  intermediate  modulus  carbon  fibers  (IM7)  and  a 
thermoplastic  polyimide  matrix  (K3B).  Unidirectional  plates  of  IM7/K3B 
composite  were  obtained  from  Boeing,  St.  Louis.  Longitudinal,  transverse  and  in¬ 
plane  shear  tests  were  conducted  under  mo  no  tonic  loading.  Typical  stress-strain 
curves  are  shown  in  Figs.  1-3.  Results  are  summarized  in  Table  1. 


Fig.  1.  Stress-Strain  Curves  of  Unidirectional  IM7/K3B 
Composite  under  Longitudinal  Tension 
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Fig.  3.  Shear  Stress-Strain  Curve  of  Unidirectional 
IM7/K3B  Composite  under  In-Plane  Shear 

Table  1  Properties  of  Unidirectional  Carbon/Polyimide  Composite 

(IM7/K3B) 


Longitudinal  Modulus,  Ei 

179  GPa  (25.9  Msi) 

Transverse  Modulus,  E2 

9.45  GPA  (1.37  Msi) 

In-Plane  Shear  Modulus,  G12 

7.93  GPa  (1.1 5  Msi) 

Major  Poisson’s  Ratio,  V12 

0.29 

Minor  Poisson’s  Ratio  V21 

0.015 

Longitudinal  Tensile  Strength,  Fit 

3,060  MPa  (444  ksi) 

Transverse  Tensile  Strength,  F2t 

71  MPa  (10.3  ksi) 

In-Plane  Shear  Strength,  F6 

88  MPa  (12.8  ksi) 

Ultimate  Longitudinal  Tensile  Strain, 

0.016 

Ultimate  Transverse  Tensile  Strain, 

0.008 

Ultimate  In-Plane  Shear  Strain, 

0.056 

Glass  Transition  Temperature,  Tg 

240°  C  (464°  F) 

Carbon/Epoxy  (EM7/977-3) 

Carbon/epoxy  prepreg  (IM7/977-3)  was  obtained  in  unidirectional  form 
and  processed  into  laminates  by  the  autoclave  process.  The  originally 
recommended  curing  cycle  by  the  manufacturer  (Cytek/Fiberite)  produced  high 
residual  stresses  which  resulted  in  low  values  for  transverse  tensile  strength  as 
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well  as  early  transverse  matrix  cracking  in  crossply  laminates.  Tests  were 
conducted  on  the  neat  resin  (977-3)  by  means  of  differential  scanning  calorimetry 
(DSC)  to  measure  the  degree  of  cure  (conversion)  as  a  function  of  time  at 
different  constant  temperatures.  Results  of  such  tests  are  shown  in  Fig.  4. 


Fig.  4.  Conversion-Time  Data  at  Different  Temperatures  for 
Cytek-Fiberite  977-3  Resin 


The  kinetic  model  for  conversion  of  the  polymer  is  given  as  [1]: 


where  P  =  conversion  fraction 
A,  E,  R,  M,  n  =  constants 
T  =  absolute  temperature 


(1) 
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The  above  was  simplified  by  fitting  an  exponential  form  to  the  DSC  data 
of  Fig.  4  as  follows: 

P  =  l-e“at  (2) 

where  a  =  a  (T)  function  of  temperature 
t  =  time 

and  the  conversion  rate  is 

—  =  ae“at  =  a(l-P)  (3) 

dt 

For  an  arbitrary  curing  cycle  T  (t)  shown  in  Fig.  5,  the  degree  of  cure  at 
timet  is  given  by 

P(t)={  A(t)e-A<’>dt=P(Te,te)  (4) 

where  A(t)  =  a  [T  (t)] 

Te,te  =  equivalent  constant  temperature  and  time  for  same  degree  of  cure 
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Residual  stresses  are  a  consequence  of  differential  thermal  shrinkage  of 
the  resin  during  cure.  Shrinkage  consists  of  two  parts,  chemical  (associated  with 
crosslinking)  and  thermal.  Chemical  shrinkage  is  irreversible,  whereas  thermal 
shrinkage  is  reversible  and  becomes  zero  at  the  “stress  free”  temperature.  The 
major  part  of  residual  stresses  is  assumed  to  be  thermal  and  thus  this  residual 
stress  is  approximately  proportional  to  the  difference  between  the  effective  stress- 
free  temperature,  Te,  and  the  reference  temperature,  T0  as 
aRoc(Te-T0) 

For  a  multi-step  curing  cycle,  where  each  step  corresponds  to  a  temperature 
increment  of  AT;  and  a  conversion  fraction  P, ,  the  thermal  residual  stress  is 
proportional  to  the  sum 

cR«£(AT),P,  (6) 

z= 1 

where  k  =  number  of  heating  steps  in  the  cycle 
k 

P,-  =  1  for  complete  cure 

;=i 

Several  two-step  curing  cycles  were  evaluated  by  fabricating 
antisymmetric  [02/902]  laminates  and  measuring  the  resulting  warpage  as 
shown  in  Fig.  6.  The  measured  and  predicted  relative  reductions  in  residual 
stresses  for  the  three  curing  cycles  tested  are  in  qualitative  agreement  as  shown  in 
Fig.  7. 
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Fig.  6.  Two-Step  Curing  Cycles  and  Contour  Maps  of 
Warped  Antisymmetric  [O2/9O2]  Laminates 
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Fig.  7.  Normalized  Residual  Stress  as  a  Function  of 
First  Curing  Step  Duration 

Unidirectional  laminates  were  fabricated  using  the  improved  curing  cycle 
(No.  3  in  Fig.  6).  The  curing  cycle  adopted  as  illustrated  in  Fig.  8  consists  of  the 
following  steps: 

1.  Place  layup  in  autoclave. 

2.  Apply  vacuum  of  at  least  85  kPa  (25  in.  Hg) 

3.  Apply  pressure  of  690  ±  35  kPa  (100  ±  5  psi)  and  release  vacuum. 

4.  Raise  temperature  from  ambient  to  280°  F  at  a  rate  of  4°  F/min. 

5.  Hold  temperature  at  280°  F  ±  10°  F  under  pressure  of  690  ±  35  kPa 

(100  ±  5  psi)  for  15  hours. 

6.  Raise  temperature  to  350°  F  ±  10°F  at  a  rate  of  4°  F/min. 

7.  Hold  temperature  at  350°  F  ±  10°  F  under  pressure  of  690  ±  35  kPa 

(100  ±  5  psi)  for  6  hours. 

8.  Cool  down  to  room  temperature  at  a  rate  of  4°  F/min.  Pressure  can  be 
released  after  the  temperature  drops  below  100°  F. 

9.  Remove  plate  from  autoclave. 
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Fig.  8  Cure  Schedule  for  IM7/977-3  Carbon/Epoxy  Composite 

Standard  characterization  tests  were  performed  on  unidirectional 
specimens  at  three  temperatures,  24°  C  (75°  F),  93°  C  (200°  F)  and  149  C  (300  F). 
Typical  stress-strain  curves  are  shown  in  Figs.  9-14.  The  properties  along  the 
fiber  direction  are  not  sensitive  to  temperature,  in  contrast  to  the  matrix- 
dominated  transverse  and  in-plane  shear  properties. 


Fig.  9.  Stress-Strain  Curve  for  [O4]  IM7/977-3  Laminate  under  Uniaxial 
Tensile  Loading  (Room  Temperature) 
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Fig.  10.  Stress-Strain  Curve  for  [908]  IM7/977-3  Laminate  under  Uniaxial 
Tensile  Loading  (Room  Temperature) 


Fig.  11.  Shear  Stress  versus  Shear  Strain  in  [106]  IM7/977-3  Specimen  under 
UniaxialTensile  Loading  (Room  Temperature) 
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Fig.  12.  Stress-Strain  Curve  for  [04]  IM7/977-3  Laminate  under  Uniaxial 

Tensile  Loading  (149°  C,  300°  F) 


Fig.  13.  Stress-Strain  Curve  for  [90s]  IM7/977-3  Laminate  under  Uniaxial 
Tensile  Loading  (149°  C,  300°  F) 
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Fig.  14.  Shear  Stress  versus  Shear  Strain  for  [106]  IM7/977-3  Specimen  under 
Uniaxial  Tensile  Loading  (149°  C,  300°  F) 

The  thermal  behavior  of  a  unidirectional  lamina  can  be  fully  characterized 
in  terms  of  two  principal  coefficients  of  thermal  expansion  cti  and  012. 
Determination  of  these  coefficients  consists  of  measuring  the  corresponding 
thermal  strains  in  a  unidirectional  composite  specimen  as  a  function  of 
temperature.  The  strain  gage  method  has  been  shown  to  be  a  practical  and 
adequate  means  of  measuring  thermal  strains  in  composites  [2].  Measured 
thermal  strains  are  plotted  versus  temperature  as  shown  in  Fig.  15.  The  slope  of 
the  curve  gives  the  coefficient  of  thermal  expansion. 

All  measured  properties  of  the  unidirectional  IM7/977-3  carbon/epoxy  composite 
are  summarized  in  Table  2. 
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Fig.  15.  Thermal  Strains  versus  Temperature  in  Unidirectional 

IM7/977-3  Composite 

Carbon/Polyimide  (IM7/RP46) 

This  material  was  obtained  in  prepreg  form  from  NASA-Langley 
Research  Center.  The  matrix  is  a  thermoset  polyimide  containing  a  solvent.  The 
recommended  curing  cycle  was  followed  [3].  Results  obtained  from  the  first 
batch  of  material  fabricated  are  summarized  in  Table  3.  The  fiber  volume  ratio 
and  other  mechanical  properties  are  relatively  low,  although  the  material  is  said  to 
have  ahigh  glass  transition  temperature  (Tg=  400°  C). 
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Table  2  Properties  of  Unidirectional  IM7/977-3  Carbon/Epoxy  Composite 


Temperature  24°C  (75°  F) 

93  °C  (200°  F) 

149°C  (75°  F) 

Longitudinal  modulus, 

Ei,  GPa(Msi) 

191  (27.7) 

186  (27.0) 

179  (26.0) 

Transverse  modulus, 

E2,  GPa  (Msi) 

9.94  (1.44) 

9.48  (1.38) 

8.93  (1.30) 

In-plane  shear  modulus, 

Gn,  GPa  (Msi) 

7.79  (1.13) 

7.45  (1.08) 

6.53  (0.95) 

Major  Poisson’s  ratio, 

Vl2 

0.35 

0.35 

0.34 

Minor  Poisson’s  ratio, 

V21 

0.018 

0.018 

0.017 

Longitudinal  tensile  strength. 
Fit,  MPa  (ksi) 

3250  (472) 

3167(460) 

2961 (430) 

Ultimate  longitudinal 

tensile  strain,  e“t 

0.0156 

0.0150 

0.0139 

Transverse  tensile  strength, 

F2t,  MPa  (ksi) 

62  (8.9) 

55  (7.9) 

51  (7.3) 

Ultimate  transverse  tensile 
strain,  s^t 

0.0060 

0.0058 

0.0056 

In-plane  shear  strength, 

F6,  MPa  (ksi) 

75  (10.9) 

61  (8.9) 

50  (7.2) 

Ultimate  in-plane  shear 

Strain,  Y6 

0.0140 

0.0132 

0.0121 

Longitudinal  thermal  expansion 
coefficient,  ai,  ps/°C  (ps/°F)  -0.9  (-0.5) 

-0.9  (-0.5) 

-0.9  (-0.5) 

Transverse  thermal  expansion 
coefficient,  <x2,  ps/°C  (ps/°F) 

22.3  (12.4) 

23.2  (12.9) 

25.9  (14.4) 

Glass  transition  temperature, 

Tg,  °C  (°F) 

150°C  (302°F) 
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Table  3  Properties  of  Unidirectional  Carbon/Polyimide  (IM7/RP46) 


Fiber  volume  ratio,  Vf 

0.53 

Longitudinal  Modulus,  Ei 

135  GPa  (19.6  Msi) 

Transverse  Modulus,  E2 

7.7.  GPa  (1.11  Msi) 

In-Plane  Shear  Modulus,  G12 

5.2  GPa  (0.75  Msi) 

Major  Poisson’s  Ratio,  V12 

0.32 

Longitudinal  Tensile  Strength,  Fit 

1,656  MPa  (240  ksi) 

Transverse  Tensile  Strength,  F2t 

52  MPa  (7.5  ksi) 

In-Plane  Shear  Strength,  F6 

83  MPa  (12.1  ksi) 

Ultimate  Longitudinal  Tensile  Strain,  e  ^ 

0.012 

Ultimate  Transverse  Tensile  Strain,  s  2t 

0.007 

Glass  Transition  Temperature,  Tg 

400°  C  (750°  F) 
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PHYSICAL  AGING  AND  CREEP 


Measurement  and  Modeling  of  Physical  Aging 

It  has  been  known  that  amorphous  materials  and  especially  polymers  are 
not  in  thermodynamic  equilibrium  below  their  glass  transition  temperature  (Tg). 
Over  a  temperature  range  below  Tg  they  tend  to  approach  the  equilibrium  state  in 
a  process  called  physical  aging.  During  this  process  the  material  becomes  more 
glasslike  with  accompanying  increase  in  density,  brittleness  and  stiffness  and  a 
decrease  in  creep  compliance.  Whereas  physical  changes  may  be  small,  the 
effects  of  aging  on  time  dependent  behavior,  such  as  creep,  can  be  very 
pronounced.  Extensive  research  has  been  reported  on  the  nature  and  effects  of 
physical  aging  of  polymers  starting  with  the  pioneering  works  of  Kovacs  [4]  and 
Struik  [5].  The  latter  determined  the  effects  of  aging  on  viscoelastic  behavior  by 
conducting  a  sequence  of  short-time  (momentary)  creep  tests  during  the  aging 
process.  He  introduced  the  concept  of  an  aging-induced  shift  factor  and  produced 
creep  master  curves.  The  effects  of  physical  aging  on  creep  of  thermoplastic 
composites  were  studied  more  recently  by  Nichols  et  al.  [6],  Sullivan  [7],  Hastie 
and  Morris  [8],  Skontorp  and  Wang  [9]  and  Gates  and  Feldman  [10]. 

Although  the  interaction  of  physical  aging  and  viscoelastic  behavior  seems 
to  be  well  understood  and  can  be  determined  under  isothermal  conditions,  it  is 
difficult  to  apply  results  to  prediction  of  long  term  behavior  of  composites 
undergoing  environmental  and  loading  fluctuations.  Given  a  material  of  unknown 
aging  history  it  would  be  desirable  to  assess  its  age  and  its  effect  on  future 
behavior.  It  would  be  helpful  if  the  “degree”  of  aging,  instead  of  just 
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chronological  age,  were  quantified  and  then  related  to  viscoelastic  effects.  The 
known  physical  change  in  a  polymer  during  physical  aging  is  a  reduction  in  free 
volume  which  results  in  volumetric  shrinkage.  The  objective  of  this  work  is  to 
measure  directly  this  aging-induced  volumetric  shrinkage  and  propose  a  model  to 
quantify  physical  aging.  Techniques  were  used  for  direct  measurement  of  aging 
strains  and  separation  of  creep  and  aging  effects  [11]. 

The  material  investigated  was  IM7/K3B,  a  carbon  fiber  reinforced 
thermoplastic  polyimide.  This  material  undergoes  physical  aging  below  its  glass 
transition  temperature  of  Tg  —  240°  C.  The  material  was  available  in  the  form  of 
24-ply,  3.15  mm  thick  unidirectional  plates.  Unidirectional  coupons  11.4  cm 
(4.5  in.)  long  and  1.27  cm  (0.5  in.)  wide  were  prepared  and  instrumented  with 
high  temperature  strain  gages  oriented  in  the  transverse  to  the  fiber  direction 
(WK-00-250BG-350,  Micromeasurements)  as  shown  in  Fig.  16. 
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Fig.  16.  Unidirectional  Composite  Specimen  with  Strain  Gage  in  the 
Transverse  (2)  Direction 
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The  volumetric  shrinkage  experienced  by  the  composite  due  to  aging  of 
the  matrix  is  directly  related  to  the  aging  strain  in  the  transverse  direction.  The 
total  volumetric  change  of  the  composite  is 


dVc 

Vc 


=  8]  +  S2  +  S3 


(7) 


where  Sj  +  s2  +  s3  =  strains  along  the  three  principal  material  axes.  Since  it  is 
not  expected  for  the  carbon  fiber  to  undergo  aging  and  since  the  2-3  plane  can  be 
assumed  to  be  a  plane  of  isotropy  for  the  material. 


^  =  s2  +  e332e2  (8) 

'c 

thus,  the  aging  strain  in  the  transverse  (2)  direction  is  one-half  the  volumetric 
change  of  the  composite.  Furthermore,  since  only  the  matrix  undergoes  aging,  the 
volumetric  change  of  the  matrix  is  related  to  the  measured  strain  as 


or 


dV, 


m  _ 


dV„ 


Vm  (1-Vf)vc  O-Vf) 

Jl-Vf)fdV. 

s2  - - - - 


e2 


vVmy 


where  Vm=  matrix  volume 

Vf  =  fiber  volume  ratio 


(9) 


In  order  to  correct  for  the  temperature  effects  on  the  gage  reading,  a 
titanium  silicate  specimen  instrumented  with  a  similar  gage  was  used  as  a 
reference  [2].  The  composite  specimen  was  heated  above  Tg  =  240°  F  to  “erase” 
any  prior  memory  and  then  quenched  down  to  the  selected  test  temperature  inside 
a  thermal  chamber.  Recording  of  strain  readings  from  the  composite  and 
reference  specimens  began  as  soon  as  temperature  equilibrium  was  reached. 
Tests  were  conducted  at  various  temperatures  in  the  range  160-190°  C 
(320-375°  F).  Typical  aging  strains  in  the  composite  in  the  transverse  direction 
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were  obtained  as  a  function  of  time  at  various  temperatures  (Fig.  17).  As  can  be 
seen,  the  magnitude  of  the  aging  strain  at  any  given  aging  time  increases  with 
temperature  at  least  in  the  time  range  investigated.  Furthermore,  all  aging  curves 
appear  to  reach  asymptotic  values  which  increase  with  temperature.  This  appears 
to  be  contrary  to  the  concept  of  free  volume  where  the  ultimate  reduction  in  free 
volume,  thereby  aging  strain,  should  be  higher  at  lower  temperatures  farther 
removed  from  Tg  [12].  However,  as  observed  by  Struik  [5],  the  slopes  of  volume 
versus  temperature  curves  at  fixed  aging  times  are  lower  than  that  of  the  as 
quenched  material  (zero  aging)  and  decrease  with  increasing  aging  time  (Fig.  18). 
This  is  attributed  to  “freezing”  of  processes  contributing  to  volume  relaxation 
below  Tg  [4]  and  is  consistent  with  the  present  experimental  results.  Similar 
trends  for  the  initial  stage  of  aging  have  been  observed  by  Kung  et  al.  [13]  who 
also  found  a  change  in  Tg  with  aging  time. 


Fig.  17.  Transverse  Aging  Strains  in  Unidirectional  Carbon/Polyimide  at 
Various  Temperatures  as  a  Function  of  Time 
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Fig.  18.  Schematic  Curves  of  Volume  versus  Temperature  Showing  Slope 
Reduction  with  Increasing  Aging  time 


The  typical  aging  strain  curve  shown  in  Fig.  19  was  modeled  by  the 
following  exponential  form 

g~£°0  —  c-(te- to)/x(T)  (10) 

80  —  8oo 

where  e =s  (T,  te )  =  current  aging  strain  at  time  te  and  temperature  T 
Sq  =  reference  strain  at  time  to 

eoo  =  8oo  (t)  =  equilibrium  aging  strain  at  temperature  T 
x(t)  =  relaxation  time  at  temperature  T 


21 


Fig.  19.  Exponential  Modeling  of  Physical  Aging  of  IM7/K3B 
Carbon/Polyimide  Composite 

This  form  was  corroborated  by  plotting  results  of  Fig.  17  on  semi-log 
scale  where  log  [(£-£00)7(80  -  £00  )3  seems  to  vary  nearly  linearly  with  aging 
time  (Fig.  20).  The  apparent  deviations  occur  at  long  aging  times  where  aging 
strain  changes  are  small  and  of  the  order  of  experimental  resolution. 

The  aging  rate  obtained  from  eq.  (4)  is 

—  do 

dte  x 

This  model  is  very  similar  to  the  one  proposed  by  Kovacs  for  isothermal 
volume  relaxation  [4].  The  main  difference  is  that  in  Kovacs’  model  the 
relaxation  time,  x,  is  also  a  function  of  aging  state,  in  this  case  aging  strain,  as 
x(e,  T).  However,  the  structural  relaxation  time  x  need  not  be  affected  by  aging  in 
the  same  way  as  the  viscoelastic  relaxation  time. 
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Fig.  20.  Normalized  Aging  Strains  versus  Time  at  Three  Different  Temperatures 


The  dependence  of  x  on  temperature  alone  was  modeled  as 


x  =  C 


/"  T  _  T,N\P 

_ _ 

V  Tg  ) 


(12) 


where  C  and  p  are  constants.  A  good  fit  of  this  form  to  experimental  results  is 
shown  in  Fig.  21.  It  was  determined  that  p  =  1.85  for  the  material  studied  here. 

Assuming  the  model  described  above  to  be  valid,  it  would  be  possible  to 
describe  aging  under  nonisothermal  conditions.  For  a  given  temperature  history 
T(te )  the  resulting  aging  can  be  calculated  from  the  rate  relation  (11) 


de  s  eQ 


+  -  =  ■ 


dte  x  x 


(13) 
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Fig.  21 .  Power  Law  Fit  for  Relaxation  Time 


where  the  relaxation  time  t[T(te)]  and  equilibrium  aging  strain  s00[T(t00)]are 
functions  of  temperature  and  thereby  functions  of  time.  The  solution  of  the 
differential  eq.  (13)  is 


1  f  fte  e« 

e  =  -  1  p— 

PL^O  T 

J*te  1 

where  p=  exp  -dte 
L^OT 


•dte  +  s0 


The  analysis  above  can  be  applied  to  the  specific  case  illustrated  in  Fig. 
22,  where  the  sample  was  first  aged  at  177°C  for  a  certain  period  of  time,  then 
cooled  down  to  149°C  and  aged  for  another  period  of  time  and  then  aged  at  177°C 
for  another  time  period.  An  isothermal  aging  curve  for  177°C  was  constructed  as 
follows.  The  segment  corresponding  to  149°C  was  transformed  on  the  time  axis 
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Fig.  22.  Construction  of  Master  Isothermal  Aging  Curve  from  Segments  of 

Isothermal  Aging  Curves 

and  shifted  to  fit  the  curve  of  177°C  aging.  The  last  segment  at  177°  C  needs 
only  to  be  shifted  along  the  time  axis. 

The  techniques  developed  give  direct  quantitative  measurement  of  aging 
strain,  thereby  volumetric  shrinkage  due  to  aging.  It  might  be  questioned  whether 
the  aging  strains  recorded  might  be  affected  or  even  overshadowed  by  other 
effects,  such  as  relaxation  of  residual  stresses.  Such  effects,  however,  would  have 
produced  positive  strains  instead  of  negative  and  would  be  higher  at  lower 
temperatures  where  the  residual  stresses  are  higher. 

The  aging  strains  measured  are  negative  and  are  higher  at  higher 
temperatures.  The  latter  trend  appears  to  be  contrary  to  the  concept  of  free 
volume,  according  to  which  the  ultimate  volumetric  reduction  should  be  higher  at 
lower  temperatures.  There  is  no  experimental  evidence  for  this  trend  and 
explanations  have  been  offered  for  the  opposite  trend,  i.e.,  the  one  observed  here 

[4]. 
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Creep  Behavior 

Creep  tests  were  conducted  on  the  same  type  of  specimen  used  for 
measurement  of  aging  [11].  The  specimens  were  24-ply  unidirectional  IM7/K3B 
coupons  instrumented  with  the  same  type  of  strain  gages  (WK-00-250BG-350) 
and  a  thermocouple.  They  were  loaded  under  a  transverse  tensile  stress  of 
approximately  30%  of  the  transverse  strength  of  the  specimen  at  the  test 
temperature.  The  creep  specimen  was  placed  in  an  environmental  chamber  along 
with  a  similarly  instrumented  dummy  specimen  (Fig.  23).  The  strain  readings 
from  the  dummy  specimen,  representing  aging,  were  subtracted  from  the  total 
readings  from  the  active  specimen  to  obtain  the  pure  creep  strain  (Fig.  24). 


Fig.  23.  Creep  Testing  Setup  in  Environmental  Chamber 
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Fig.  24.  Isothermal  Creep  Test  of  IM7/K3B  at  Elevated  Temperature  (350°  F) 

The  effect  of  aging  on  creep  was  checked  at  one  temperature  only,  149°  C, 
by  means  of  short  time  creep  tests  at  different  aging  times.  No  noticeable  effects 
were  detected,  primarily  due  to  the  small  amount  of  creep  and  small  degree  of 
aging.  Despite  the  low  creep  recorded,  it  was  noticed  that  the  material 
experienced  some  degradation  with  time.  The  streaks  shown  in  Fig.  25  did  not 
cause  failure  but  appear  to  be  precursors  of  transverse  matrix  cracks. 


Fig.  25.  Damage  Mechanisms  during  Creep  of  IM7/K3B  under 
Transverse  Tension 
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Nonisothermal  creep  behavior  was  briefly  investigated  in  a  qualitative 
manner.  The  average  creep  compliance  was  determined  under  cyclic  thermal 
conditions  with  the  temperature  cycled  between  163  and  190°  C  and  compared 
with  isothermal  creep  compliances  at  177  and  190°  C  (Fig.  26).  Creep  under 
cyclic  thermal  conditions  is  higher  than  isothermal  creep  at  the  mean  cyclic 
temperature  (177°  C).  This  phenomenon  could  become  more  pronounced  at 
higher  temperatures  and  under  higher  loads. 


Fig.  26.  Comparison  of  Isothermal  and  Thermal  Cycling 
Creep  Tests  of  IM7/K3B 

Pulse  Loading  Method  for  Viscoelastic  Characterization  of  Aging  Effects 

As  mentioned  before,  physical  aging  in  a  carbon/polimide  composite  was 
measured  and  modeled.  Whereas  physical  aging  changes  may  be  small,  the 
effect  of  aging  on  time-dependent  behavior,  such  as  creep,  can  be  very 
pronounced.  This  is  usually  evaluated  by  conducting  a  sequence  of  short-time 
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(momentary)  creep  tests  during  an  isothermal  aging  process.  Such  tests  are  time 
consuming  and  may  be  affected  by  the  aging  process  itself. 

A  dynamic  pulse  loading  method  was  proposed  for  characterizing  the 
viscoelastic  behavior  of  the  material  at  any  instant  during  its  aging  process  [14]. 
The  method  is  based  on  dynamic  pulse  loading  of  short  duration  which  implies 
quick  recovery.  It  is  best  suited  for  tensile  loading  and  yields  convenient 
viscoelastic  parameters,  i.e.,  complex  modulus.  It  is  inexpensive  and  simple  to 
implement.  An  instrumented  drop  tower  apparatus  was  designed  and  built  to 
apply  a  dynamic  pulse  loading  to  a  specimen  (Fig.  27).  The  specimen  loading 
pulse,  measured  with  a  load  cell,  and  the  specimen  strain  pulse  measured  directly 
with  strain  gages  are  analyzed  to  yield  the  complex  modulus  as  a  function  of 
frequency.  The  instrumentation  set  up  is  shown  in  Fig.  28. 

A  sequence  of  pulse  loadings  was  applied  during  the  aging  process  as 
shown  schematically  in  Fig.  29.  For  each  pulse  loading,  stresses  c(t)  and  strains 
s(t)  were  acquired  digitally  and  their  Fourier  transforms  were  determined.  The 
complex  modulus  E*  (z'ca)  was  calculated  as 

E*(to)=E'(«,)+iE"(»)=||||  (15) 

where  the  storage  modulus  E7  (co)  and  the  loss  modulus  E"(co)  were  determined 
and  used  to  characterize  the  viscoelastic  behavior  of  the  material  over  the 
frequency  range  covered  by  the  pulse  loading. 

Each  pulse  loading  yields  results  over  a  limited  frequency  range.  The 
method  is  more  reliable  at  low  frequencies,  less  than  100  Hz.  To  improve  the 
accuracy,  multiple  hits  are  necessary.  The  method  was  applied  to  the  IM7/K3B 
carbon/polyimide  composite  undergoing  aging  under  isothermal  and 
nonisothermal  conditions  (conditions  A  and  B  in  Fig.  30). 
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Fig.  27.  Experimental  Setup  of  Pulse  Loading  Method  for  Viscoelastic  Characterization  during  Physical  Aging  Process 


Fig.  28.  Instrumentation  Setup  for  Pulse  Loading  Method 


Results  of  storage  and  loss  moduli  as  a  function  of  aging  time  for  the 
isothermal  case  A  are  superimposed  on  measured  results  of  aging  strain  for  an 
IM7/K3B  carbon/polyimide  composite  in  Fig.  31.  It  is  seen  that  there  is  a  direct 
correspondence  between  measured  aging  strain  (volumetric  shrinkage)  and 
viscoelastic  parameters.  It  is  proposed  that  volumetric  shrinkage,  as  measured 
directly  and  modeled,  is  a  true  indicator  of  the  degree  of  aging  and  not  just 
chronological  aging  and  it  can  be  used  under  both  isothermal  and  nonisothermal 
conditions.  With  a  single  measurable  and  possible  to  model  parameter 
representing  the  degree  of  physical  aging,  the  next  question  is  whether  the 
viscoelastic  behavior  is  uniquely  related  to  the  degree  of  aging  regardless  of  the 
prior  thermal  history. 
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Storage  Modulus,  E’(a>),  (MPa) 


In  the  case  of  nonisothermal  aging,  results  of  aging  strain  and 
corresponding  storage  and  loss  modulus  can  be  converted  by  time  shifting  into 
results  for  any  selected  isothermal  reference  condition.  Thus,  results  for 
temperature  history  (B)  of  Fig.  30,  can  be  converted  to  equivalent  results  for  the 
isothermal  temperature  history  (A)  as  shown  in  Fig.  32. 

The  proposed  pulse  loading  method  is  intended  to  replace  the  sequence  of 
short  time  creep  tests  by  a  sequence  of  pulse  loading  tests  during  the  aging 
process.  The  purpose  was  to  provide  a  faster,  more  representative  and  convenient 
way  to  evaluate  the  viscoelastic  behavior  for  materials  during  the  aging  process. 
By  observing  the  change  of  storage  and  loss  moduli  it  is  easier  for  one  to  monitor 
the  change  of  viscoelastic  behavior  of  materials  during  the  aging  process,  and 
hopefully  shed  new  light  into  the  relationship  between  aging  and  viscoelastic 
behavior.  The  results  obtained  to  date  show  that  there  is  a  direct  correlation 
between  viscoelastic  properties  and  degree  of  aging.  This  implies  that  there 
should  be  a  way  to  assess  the  current  age  of  a  material  with  unknown  thermal 
history  and  predict  its  long  term  behavior. 
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Fig.  32.  Conversion  of  Nonisothermal  Aging  Results  to  Equivalent  Isothermal  Results 


TRANSVERSE  CRACKING  IN  CROSSPLY  LAMINATES 
UNDER  MONOTONIC  TENSION 

Introduction 

Damage  in  multidirectional  laminates  consists  of  the  fundamental  failure 
modes  of  matrix  cracking  in  the  off-axis  plies,  delaminations  and  fiber  failures 
[15-18].  Experimental  observations  have  shown  that  the  first  stage  of  damage  is 
transverse  matrix  cracking  in  the  90°  plies  [15,  19,  20].  Although  the  formation 
of  transverse  cracks  by  itself  does  not  cause  catastrophic  failure,  its  presence  can 
influence  the  overall  mechanical  behavior  of  the  laminate,  such  as  reduction  of 
the  laminate  stiffiiess.  Moreover,  the  existence  of  transverse  cracks  degrades  the 
strength  of  composite  laminates  by  inducing  other  more  severe  damage  modes. 

Two  main  approaches  have  been  proposed  to  analyze  the  transverse 
cracking  problem.  One  of  these  is  based  on  lamina  strength  [21, 22,  30,  31]  while 
the  other  is  based  on  fracture  mechanics  [23,  24,  31].  In  the  fracture  mechanics 
approach  a  microcrack  will  form  when  the  strain  energy  release  rate  reaches  some 
critical  value.  However,  this  approach  requires  the  determination  of  the  critical 
fracture  toughness,  Gc  ,  with  additional  not  easy  to  conduct  experiments. 
Moreover,  out  experiments  show  that  for  the  IM7/977-3  composite  laminate 
almost  all  transverse  cracks  propagate  through  the  width  of  the  transverse  plies 
nearly  instantaneously  during  the  test,  which  means  that  the  lamina  strength 
model  is  precise  enough  when  dealing  with  this  problem. 

Reifsnider  et  al.  [25]  predicted  the  stiffiiess  reduction  using  a  shear  lag 
analysis.  Dvorak  and  Laws  [26]  proposed  a  progressive  damage  analysis  based 
on  statistical  fracture  mechanics  for  transverse  cracking  in  crossply  composites. 
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This  model  is  based  on  a  shear  lag  model  which  gives  a  good  prediction  for 
stiffness  reduction.  In  order  to  analyze  progressive  cracking,  they  introduced  the 
probability  density  function  for  the  site  of  the  next  crack.  However,  they  did  not 
proceed  to  study  the  strength  distribution  of  the  transverse  layers. 

Berthelot  and  Le  Corre  [30]  presented  a  statistical  analysis  of  the 
transverse  cracking  and  delamination  process  induced  in  crossply  laminates 
subjected  to  tensile  loading.  They  used  a  generalized  stress  model  and  statistical 
distribution  of  strength  in  the  90°  layer  to  predict  transverse  cracking.  However, 
they  only  studied  damage  at  room  temperature. 

Lee  and  Daniel  [22]  employed  a  modified  shear  lag  analysis  to  develop  a 
progressive  damage  scheme  for  crossply  composite  laminates.  The  shear  lag 
parameter  is  expressed  only  as  a  function  of  shear  moduli  of  the  lamina  and 
geometry  of  the  laminate.  They  gave  closed  form  solutions  for  the  transverse 
crack  density,  stress  distributions  and  reduced  stiffness  of  damaged  plies  as  well 
as  the  entire  laminate  as  a  function  of  applied  load  and  lamina  properties.  For 
progressive  transverse  cracking,  the  approach  assumes  the  next  set  of  cracks 
occurs  when  the  maximum  axial  stress  in  those  plies  reaches  the  transverse 
strength  of  the  layer.  The  theory  assumes  uniform  transverse  strength  along  the 
specimen,  and  in  this  case  the  new  cracks  always  occur  in  the  middle  of  the 
segments.  However,  defects  always  exist  in  the  material,  which  degrade  its 
strength.  Experiments  show  that  the  transverse  strength  of  a  unidirectional  lamina 
is  not  uniform,  but  follows  some  statistical  distribution.  Thus,  it  is  more  realistic 
to  study  progressive  transverse  cracking  considering  statistical  strength 
distribution. 
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The  shear  lag  model  proposed  by  Lee  and  Daniel  [22]  was  applied  to 
IM7/977-3  carbon/epoxy  considering  the  transverse  strength  distribution  obtained 
from  experiments.  The  progressive  damage  development  of  the  crossply  laminate 
was  studied  systematically  by  the  stochastic  method.  Monte  Carlo  simulation  was 
used  to  predict  transverse  cracking  at  both  room  and  high  temperatures. 
Transverse  cracks  were  recorded  by  X-ray  radiography  and  in-test  photography. 
A  master  curve  relating  the  normalized  transverse  stress  with  normalized  crack 
density  for  different  lay-ups  and  temperatures  was  established. 


Transverse  Cracking  -  Deterministic  Model 

In  the  previously  developed  shear  lag  model,  Lee  and  Daniel  gave  closed 
form  expressions  for  local  stress  distributions,  crack  density  and  reduced  stiffness 
in  a  crossply  laminate  as  a  function  of  applied  stress  [22].  The  axial  (x-axis) 
stress  distributions  in  the  0°  and  90°  layers  of  a  laminate  segment  containing  two 


cracks  are  given  as  (Fig.  33). 

z 


Fig.  33.  Element  of  Crossply  Laminate  with  Transverse  Cracks  in  90°  Layer 
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where 

Ei,  E2  =  Longitudinal  and  transverse  lamina  moduli 

E„  =  +^2^2  laujjnate  modulus  (before  damage) 

x  h,h2E1E2 

G12,  G23  =  in-plane  and  out-of-plane  lamina  shear  moduli 

Cx’CT2x>CT2r  =  applied  stress,  90°  layer  stress,  and  90°  layer  residual 
stress,  respectively 

a2  _  3(ht  +  h2)Ex  *  G12G23 

h1h2E1E2  h,G23+h2G12 


The  stiffness  reduction  caused  by  the  matrix  cracking  was  calculated 
based  on  the  stress  distributions  above  [22].  The  modulus  reduction  ratio  is 
expressed  as  a  function  of  material  and  geometric  parameters  and  the  crack 


spacmg  as 


,  2  E2  h2  A  ,  al  alr  Ex 

1+ - Manh —  +  -JL  — i 

al  Ejh,  2  ax  E! 


—  tanh— 
al  2 


where  Ex ,  Ex  =  initial  and  reduced  moduli  of  the  laminate,  respectively. 

Crack  density  can  be  predicted  in  a  deterministic  fashion  by  assuming 
uniform  crack  spacing  and  equating  the  value  of  a2x  in  eq.  (17)  at  x=//2  to  the 
transverse  lamina  strength  F2t 
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MX  =  -  =F2t 


At  a  value  of  a2x  slightly  less  than  F2t  the  crack  spacing  is 


,12  _i , 

I  =  —  =  —  COS  4 

X,  a 


°2x  =  F2t  -  8 


where  8  «  F2t 


c  _  ct2x 
°  2x  -  F2t 


F 

ct2x  =  =2-  ctx  +  <r2r=  average  stress  in  90°  layer  segment  between  cracks 

Ex 

At  a  value  of  cr2x  slightly  higher  than  F2t  a  new  crack  forms  in  the  middle 
between  the  existing  cracks  and  the  new  crack  spacing  is  half  the  previous  one. 

/' =  ^7  =  -  cos-1  ^  ,  or2x=F2t  +  5  (21) 

X  a 


The  most  probable  value  of  crack  density  at  a  stress  value  of  a2x  =  F2t  is 
the  value  corresponding  to  the  average  of  the  two  values  given  by  eqs.  (20)  and 


2\l  l )  2  4cosh  1  % 


or,  solving  for 


£,  =  cosh- 
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and  solving  for 


°2x 

F2t 


(24) 


(25) 


where  cr2x  =  a2x  /  F2t  and  A,  =  — ^ 

Equation  (25)  above  is  a  generalized  relation  between  normalized  applied 
stress  and  normalized  average  crack  density  and  is  independent  of  crossply  layup 
and  transverse  lamina  strength  (temperature). 


Transverse  Cracking  -  Statistical  Model 

Equation  (17)  gives  the  axial  stress  distribution  in  the  transverse  layer. 
When  the  applied  stress  is  low,  which  means  that  the  crack  spacing  is  much  larger 
than  the  thickness  of  the  transverse  layer,  the  axial  stress  in  the  central  region  of 
the  transverse  layer  is  nearly  uniform.  However,  because  of  the  existence  of 
flaws  in  the  transverse  layer,  the  transverse  strength  between  two  existing  cracks 
is  not  unique  but  falls  within  a  band  as  shown  in  Fig.  34.  With  increase  of  the 
applied  load,  the  point  where  the  stress  first  reaches  its  corresponding  strength 
value  in  the  band  will  be  the  location  of  the  next  crack.  Since  the  transverse 
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strength  varies  from  point  to  point,  the  new  crack  could  fall  anywhere  except  in 
the  two  end  regions,  the  location  of  the  crack  being  dependent  on  the  stress  and 
strength  distributions.  This  means  that  the  new  crack  does  not  necessarily  occur 
at  the  midpoint.  A  relatively  minor  statistical  variation  in  transverse  strength 
distribution  will  lead  to  a  statistical  variation  in  crack  spacing.  Consequently,  the 
transverse  strength  distribution  is  needed  to  further  investigate  transverse  cracking 
in  the  laminate. 


Fig.  34  Stress  Distribution  and  Strength  Range  Band  between  two  Existing 

Cracks  in  90°  Layer 


The  Weibull  distribution  function  was  used  to  describe  the  transverse 
strength  distribution.  The  probability  of  failure  of  a  uniformly  loaded  specimen 
predicted  by  the  Weibull  distribution  is 


P(cr)  =  1  -  exp 


V0 


f 

cr-a„  ' 


<*o  ) 


(26) 
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where  m  is  the  shape  parameter,  cy0  is  a  scale  parameter,  au  is  the  strength 
threshold,  V  is  the  specimen  volume,  and  Vo  is  reference  volume.  The  parameters 
m,  co  and  cu  can  be  determined  by  experiment. 

A  Monte  Carlo  simulation  was  used  to  analyze  transverse  cracking  in 
crossply  laminates  based  on  the  statistical  (Weibull)  distribution  of  transverse 
lamina  strength.  Given  the  strength  distribution  and  a  criterion  for  cracking,  a 
computer  simulation  is  performed  to  analyze  the  progressive  cracking  in  the 
transverse  ply. 

The  simulation  performed  in  this  study  differed  from  those  of  other 
studies,  such  as  Wang’s  [27].  Instead  of  assuming  a  flaw  size  distribution  in  the 
transverse  layer,  we  simulated  the  crack  multiplication  by  partitioning  the  initial 
gage  length  into  5000  equal  elements  and  assigning  a  strength  value  to  each 
element  in  accordance  with  our  experimental  distribution  function.  The  new 
crack  occurs  when  the  stress  in  an  element  reaches  its  strength. 

The  gage  length  in  the  simulation  was  15.2  cm  (6  in.),  which  was  the  same 
as  in  the  tests.  The  number  of  elements  should  be  large  enough  to  ensure  a 
smooth  transverse  strength  distribution  in  the  gage  length. 

Since  the  strength  of  the  transverse  layer  has  a  statistical  distribution,  the 
new  crack  location  is  not  necessarily  in  the  middle  of  the  segment,  which  means 
that  the  crack  spacing  is  not  uniform.  Instead,  the  distribution  of  crack  spacings  is 
determined  by  the  load  level  and  transverse  strength  distribution  [28].  Cho, 
Holmes  and  Barber  [29]  studied  the  matrix  crack  spacing  problem  in  a  uniaxial 
ceramic  composite.  They  also  pointed  out  that  if  the  matrix  is  completely 


44 


homogeneous  in  strength,  the  cracks  will  always  form  at  the  midpoint  of  the 
segments.  In  this  case,  they  can  determine  the  probability  density  function  (PDF) 
of  crack  spacing,  which  is 

f(/,/0)  =  /0//2  (27) 

where  /0  is  the  crack  spacing  corresponding  to  a  specific  applied  stress  ax . 

The  PDF  of  the  transverse  strength  F2t,  assuming  a  Weibull  distribution 
with  parameter  <ru  =  0,  is  expressed  as 


The  PDF  of  lo  can  be  obtained  by  statistical  theory  and  eqs.  (17)  and  (19) 

f2('o)  =  f,[F2,('o)]^!-  C®) 

d/0 

The  PDF  of  crack  spacing  considering  the  transverse  strength  distribution  is 

f,(/)=  ff  (/,/«).  f2(/„)d/0  (30) 

Thus,  the  distribution  of  crack  spacings  can  be  obtained  in  terms  of  the 
distribution  of  transverse  strengths.  Equation  (30)  is  an  integral  equation  which 
can  be  solved  numerically  to  obtain  the  crack  spacing  distribution  f3  (/) . 

Car  Parking  Analog 

An  analogy  can  be  drawn  between  transverse  cracking  in  a  crossply 
laminate  and  car  parking  along  a  straight  line.  If  all  cars  are  of  length  a  and  the 
maximum  number  of  cars  are  parked,  the  clear  spacing  between  cars  can  vary 
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between  0  and  a,  corresponding  to  car  center  spacings  of  a  and  2a  or  car  densities 
of  1/a  and  l/2a  (Fig.  35).  Kimber  and  Keer  have  solved  the  relevant  integral 
equation  of  the  car  parking  problem  and  found  that  for  a  line  of  infinite  length,  the 
average  spacing  between  car  centers  is  4a/3,  where  a  is  the  length  of  a  car  [32]. 
To  find  the  average  crack  spacing  for  multiple  matrix  cracking,  an  analogy  can  be 
made  with  the  average  spacing  between  cars  of  length  a  parked  randomly  along  a 
line.  From  this  analogy  it  follows  that  the  average  spacing  between  two  cracks  is 
4a/3,  where  a  is  the  distance  from  a  crack  to  a  point  where  the  stress  reaches  F2t. 
In  this  case,  the  average  crack  density  is  Xm  =  3/4 a .  The  distance  a  can  be 
calculated  from  the  shear  lag  model,  eq.  (17).  If  a  =  112  and  CT2X=  F2t  in  eq.  (17), 
the  same  result  is  obtained  as  in  the  deterministic  model  by  replacing  l  with  1/A,m 
in  eq.  (17). 

Experimental  Procedure 

Material  Characterization 

The  material  used  in  this  investigation  was  IM7/977-3  carbon/epoxy  and 
supplied  in  prepreg  form.  Laminated  panels  were  fabricated  using  the  autoclave 
process.  Two  unidirectional  laminates  were  made,  a  [908]  plate  for  determination 
of  transverse  properties,  and  a  [04]  plate  for  determination  of  longitudinal  and 
shear  properties.  Two  crossply  panels  were  made,  [0/904]s  and  [0/902]s.  All 
panels  were  inspected  by  ultrasonic  C-scanning  and  found  to  have  no  discernible 
defects.  The  fiber  volume  ratio  was  0.65,  which  was  determined  by  image 
analysis  of  photomicrographs  of  transverse  (to  the  fibers)  cross  sections  of  the 
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Car  parking  problem 

a  =  car  length 
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Fig.  35.  Car  Parking  Analogy  of  Transverse  Cracking  in  Crossply  Laminate 


composite  [2].  The  [90g]  panels  were  cut  with  a  diamond  saw  into  22.9  cm  (9.0 
in.)  x  2.54  cm  (1.0  in.)  coupons  for  transverse  tensile  testing.  The  [04]  panels 
were  cut  into  22.9  cm  (9.0  in.)  x  1.27  cm  cm  (0.5  in.)  coupons  with  the  fibers 
oriented  parallel  to  and  at  10°  with  the  loading  axis  for  longitudinal  tension  and 
in-plane  shear  testing,  respectively. 

To  prevent  fiber  breakage  during  gripping,  end  tabs  were  used.  For  room 
temperature  testing,  glass/epoxy  tabs  were  used.  They  were  bonded  to  the 
laminate  with  epoxy  (Hysol  907,  blue  epoxy-patch  kit,  Dexter  Corporation). 
Before  the  test,  the  epoxy  adhesive  was  cured  at  room  temperature  for  24  hours. 
For  high  temperature  tests,  aluminum  tabs  were  bonded  to  the  laminate  with 
FM300  adhesive  film  (Cytec  Fiberite  Inc.).  All  tests  were  performed  in  an  Instron 
1331  servohydraulic  testing  machine  with  an  Instron  8500  controller.  An 
environmental  chamber  attached  to  the  testing  machine  was  used  for  high 
temperature  tests. 

The  unidirectional  laminate  was  fully  characterized  at  three  different 
temperatures,  24°  C  (75°  F),  93°  C  (200°  F)  and  149°  C  (300°  F).  The  properties 
are  tabulated  in  Table  2. 

Transverse  Cracking  of  Crossply  Laminate 

Damage  development  was  studied  in  two  crossply  layups,  [0/904]s  and 
[0/902]s.  For  each  layup,  experimental  data  of  applied  stress  versus  crack  density 
were  obtained  both  at  room  temperature  24°  C  (75°  F)  and  at  149°  C  (300°  F).  In 
the  room  temperature  tests,  transverse  cracking  was  monitored  by  X-radiography 
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enhanced  with  a  penetrant  opaque  to  X-rays.  At  each  load  level,  the  specimen 
was  taken  off  the  machine  and  exposed  to  soft  X-rays  to  obtain  radiographs. 
Figure  36  shows  a  typical  progression  of  cracking  in  the  90°  layer  of  a  [0/904]s 
laminate  by  X-radiography  at  various  applied  stress  levels.  However,  in  the  high 
temperature  test  this  technique  could  not  work,  because  applying  X-rays  would 
change  the  thermal  environment  of  the  test.  Instead,  a  high  temperature  in-test 
photography  setup  was  used  to  inspect  the  damage  in  the  transverse  layer.  The 
edge  of  the  specimen  was  carefully  polished  before  the  test.  As  shown  in  Fig.  37, 
a  long  distance  microscope  (Infinity  Co.)  was  used  to  monitor  the  crack 
multiplication  from  the  edge  during  the  test.  A  Nikon  camera  connected  to  the 
microscope  was  used  to  record  the  crack  distribution  along  the  gage  length.  Two 


typical  crack  images  of  the  edge  are  shown  in  Fig.  38. 


Stress 

MPa(ksi) 


Fig.  36.  X-Radiographs  of  IM7/977-3  [0/904]s  Laminate  under  Uniaxial  Tensile 
Loading  at  Various  Applied  Stress  Levels 
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Long  distance  microscope 


Fig.  37.  Setup  for  Transverse  Crack  Measurement  at  High  Temperatures 


Fig.  38  Typical  Crack  Images  at  the  Edge  during  High 
Temperature  Test  (149°  C,  300°  F) 


Results  and  Discussion 


Statistical  Distribution  of  Transverse  Strength 

As  discussed  before,  the  transverse  strength  distribution  is  important  in 
investigating  the  progressive  damage  in  the  transverse  layer  of  the  crossply 
laminate.  Experimental  results  of  the  transverse  tensile  tests  were  plotted  in  the 
forms  of  distribution  curves  for  the  three  test  temperatures  (Fig.  39).  These 
results  were  normalized  by  dividing  the  failure  stress  by  the  corresponding  mean 
strength  (F2t)m  for  each  temperature.  It  is  seen  in  Fig.  40  that  the  normalized 
strength  distribution  is  independent  of  temperature.  Therefore,  transverse  strength 
can  be  described  by  one  set  of  Weibull  parameters,  independent  of  temperature. 
In  this  case  a  shape  parameter  m  =  16  was  found  for  all  three  temperatures. 


Fig.  39  Cumulative  Transverse  Strength  Distributions  at 
Three  Temperatures 
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Fig.  40  Cumulative  Normalized  Transverse  Strength  Distribution  at 
Three  Different  Temperatures 

Stiffness  Reduction 

The  measured  stress-strain  behavior  of  the  [0/904]s  crossply  laminate 
under  monotonic  loading  is  shown  in  Fig.  41  with  a  clear  indication  of  the  region 
of  transverse  matrix  cracking.  Similar  stress-strain  curves  under  repeated 
loading/unloading  conditions  are  shown  in  Fig.  42  where  it  is  seen  that  the 
process  of  unloading  and  reloading  does  not  introduce  any  additional  damage 
over  that  produced  by  the  monotonic  loading. 

The  masured  laminate  stiffness  as  a  function  of  transverse  crack  density 
was  compared  with  the  closed-form  expression  for  reduced  modulus  given  by  eq. 
(18)  in  Fig.  43.  The  agreement  is  very  good. 
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Stress. 


Strain,  e,  (%) 


Fig.  41.  Longitudinal  Stress-Strain  Curve  of  [(0/904)]s  IM7/977-3  Laminate 


Strain,  s,  (%) 


Fig.  42.  Stress-Strain  Curves  for  [0/904]s  Laminate  under  Monotonic 
Loading  and  Loading-Unloading  Cycles 
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Fig.  43.  Normalized  Axial  Modulus  of  [0/904]s  IM7/977-3  Laminate  as  a 
Function  of  Transverse  Crack  Density 


Transverse  Cracking 

Many  models  have  been  proposed  for  prediction  of  transverse  cracking  in 
crossply  laminates  [15,  16,  19,  22-24].  Few,  however,  discuss  the  effect  of 
laminate  layup  and  temperature.  The  generalized  expression  given  here  in  eq. 
(25)  yields  a  master  curve  relating  normalized  applied  stress  in  the  transverse 
layer  and  normalized  crack  density  independently  of  layup  and  temperature.  This 
closed  form  relation  is  in  excellent  agreement  with  experimental  results  as  shown 
in  Fig.  44.  Equation  (25)  is  a  deterministic  relation  between  loading  and  average 
crack  spacing. 
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Fig.  44.  Master  Curve  for  Transverse  Cracking  Compared  with 
Experimental  Results 

The  nonuniformity  of  crack  spacing  was  investigated  by  statistical 
analysis  taking  into  account  the  statistical  distribution  of  transverse  strengths  as 
discussed  before.  To  verify  the  validity  of  the  numerical  simulations,  results  of 
applied  stress  versus  crack  density  obtained  by  numerical  simulation  for  two 
laminates,  [0/902]s  and  [0/904]s ,  and  at  two  temperatures,  24°  C  and  149°  C,  were 
compared  with  experimental  results  in  Figs.  45  and  46.  The  agreement  is  good. 


Fig.  45.  Simulations  and  Experimental  Data  of  Applied  Laminate  Stress  versus 
Crack  Density  for  [0/904]s  Laminate  at  24°  C  and  149°  C 


55 


Fig.  46.  Simulations  and  Experimental  Data  of  Applied  Laminate  Stress  versus 
Crack  Density  for  [0/902]s  Laminate  at  24°  C  and  149°  C 

Simulations  yielded  also  crack  spacing  distributions,  which  were 
compared  with  experimental  measurements.  The  experimentally  measured  and 
predicted  by  simulation  crack  spacing  distributions  for  a  [0/904]s  laminate  under 
an  applied  stress  of  ctx  =  276MPa(40ksi)  are  shown  in  Figs.  47  and  48, 
respectively.  The  agreement  between  them  is  satisfactory. 


Crack  spacing,  £ ,  (cm) 


Fig.  47  Histogram  of  Crack  Spacings  from  Experiments 
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Crack  spacing,  £  ,  (cm) 


Fig.  48.  Histogram  of  Crack  Spacings  from  Simulation 

The  statistical  models  proposed  by  eqs.  (27)  and  (30)  were  used  to 
compare  the  probability  density  functions  of  crack  spacings  with  that  of  numerical 
simulation  as  shown  in  Fig.  49.  The  new  model  given  by  eq.  (30)  which 
incorporates  the  transverse  strength  distribution  obtained  experimentally,  agrees 
well  with  the  simulations  and  experiments.  The  approximate  model  given  by  eq. 
(27),  which  assumes  uniform  strength,  does  not  agree  well  with  the  simulations. 


Fig.  49.  Comparison  of  Simulation  and  Theoretical  Models  Regarding  PDF  of 

_  E2  _ 

Crack  Spacings  (for  ax2  =  =^-crx  +  ar2  =  1-1  ^t) 

Ex 
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The  statistical  models  using  simulations  and  the  car  parking  analogue, 
when  expressed  in  terms  of  average  crack  density,  yield  the  same  relationship 
between  applied  stress  and  crack  density  as  the  deterministic  model  discussed 
before  (eq.  25,  Fig.  44). 

Summary  and  Conclusions 

Transverse  cracking  in  a  monotonically  loaded  crossply  laminate  was 
studied  following  a  Monte-Carlo  simulation  based  on  the  statistical  distribution  of 
transverse  lamina  strengths.  It  was  found  experimentally  that  the  normalized  (by 
the  mean  value)  transverse  strength  can  be  described  by  a  Weibull  distribution 
function  with  parameters  independent  of  temperature. 

The  statistical  distribution  of  transverse  strengths  was  used  to  analyze  the 
transverse  crack  multiplication  including  a  distribution  of  nonuniform  crack 
spacings.  Simulation  results  were  in  agreement  with  experimental  ones. 

In  addition,  a  general  relationship  was  developed  between  normalized 
applied  stress  and  normalized  crack  density  based  on  a  previously  developed 
shear  lag  model  (Lee  and  Daniel  [22]).  This  relationship  is  independent  of 
crossply  laminate  layup  and  temperature.  It  was  shown  that  all  predictions  agree 
with  each  other  and  with  experimental  results. 
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DAMAGE  DEVELOPMENT  UNDER  CYCLIC  LOADING 
Introduction 

Fatigue  behavior  of  composite  materials  is  a  complex  phenomenon 
consisting  of  cyclic  load-induced  events  and  processes  which  determine  long-term 
performance  of  composites.  Thus,  a  great  deal  of  research  has  been  focused  on 
the  fatigue  damage  accumulation  and  fatigue  life  prediction  of  composite 
structures.  Although  many  theoretical  and  experimental  investigations  have  been 
reported,  providing  a  fairly  comprehensive  understanding  of  damage  mechanisms 
developed  during  fatigue  loading  of  fiber  reinforced  composite  materials,  it  is  still 
an  open  question  for  the  integrity  analysis  of  composites  to  predict  the  fatigue 
damage  and  long-term  behavior  of  a  laminate  subjected  to  specified  load  and 
environment  conditions  (temperature,  moisture,  etc.)  [33-46].  Specifically, 
temperature  is  a  very  important  factor  which  influences  fatigue  behavior  of 
composite  materials.  Many  experiments  show  that  fatigue  damage  at  elevated 
temperatures  develops  much  faster  than  at  room  temperature  [41,  46].  Thanks  to 
the  time-temperature  superposition  theory  for  viscoelastic  analysis,  the  possibility 
of  extending  the  time  scale  (number  of  cycles)  by  varying  the  temperature  is  of 
particular  interest  in  fatigue  studies. 

Broutman  and  Sahu  first  assumed  a  monotonically  decreasing  residual 
strength  with  cyclic  loading  [33].  A  modified  Miner  s  rule  was  presented  based 
on  linear  strength  reduction  as  a  function  of  fractional  life  spent  at  a  given  stress 
level. 
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Reifsnider  and  Stinchcomb  [34]  proposed  the  critical  element  concept 
which  is  capable  of  describing  and  predicting  the  strength  and  life  of  composite 
laminates  subjected  to  general  cyclic  loading.  There  are  two  main  ideas  in  this 
methodology.  First,  the  structure  is  separated  into  critical  and  subcritical 
elements.  Secondly,  the  damage  evolution  is  characterized  by  its  remaining 
strength.  The  critical  element  is  that  portion  of  the  structure  which  controls 
failure.  The  subcritical  element  does  not  result  in  structural  failure,  but  it  does 
affect  the  stress  state  in  the  critical  element.  The  model  provides  a  general 
framework  for  determining  fatigue  response  of  composite  materials  from  an 
understanding  of  microstructural  damage  and  failure  mechanisms  of  materials. 
However,  the  application  of  this  model  requires  micromechanical  formulation  of 
damage  and  damage  evolution  to  represent  the  local  load  share  and  failure  which 

determines  the  final  failure  of  the  laminate. 

Lee,  Daniel  and  Yaniv  [35]  proposed  a  deterministic  model  to  predict  the 
fatigue  life  of  crossply  laminates  based  on  the  static  and  fatigue  properties  of  the 
unidirectional  lamina.  In  their  model,  the  fatigue  life  consists  of  two  parts.  The 
first  part  corresponds  to  the  first  stage  of  damage  development  consisting  of 
transverse  matrix  cracking  and  multiplication  up  to  the  limiting  characteristic 
damage  state  (CDS).  This  part  is  determined  from  the  S-N  curve  of  90°  plies. 
The  second  part  following  attainment  of  CDS  is  based  on  the  state  of  stress  in  the 
0°  plies  and  the  S-N  curve  of  the  0°  lamina.  It  is  assumed  that  the  0°  plies  in  the 
damaged  laminate  behave  under  cyclic  loading  like  the  0°  unidirectional  lamina. 
The  total  life  predicted  is  taken  as  the  sum  of  the  two  parts. 


60 


Diao,  Ye  and  Mai  [39]  proposed  a  statistical  model  of  residual  strength 
and  fatigue  life  of  composite  laminates  on  the  basis  of  the  statistical  nature  of 
fatigue  damage  and  the  concept  of  a  representative  volume  for  damage  evolution. 
In  this  model,  the  shear  lag  model  proposed  by  Lee  and  Daniel  [22]  was  used  to 
obtain  the  stress  redistribution  function,  which  characterizes  the  stress 
redistribution  process  between  the  critical  and  subcritical  elements.  However,  the 
present  stress  redistribution  function  is  only  based  on  the  transverse  matrix 
cracking  in  the  90°  plies  up  to  the  CDS.  Other  damage  mechanisms  need  to  be 
considered  to  obtain  a  more  precise  stress  redistribution  function. 

Rotem  and  Nelson  [41]  studied  fatigue  behavior  of  graphite-epoxy 
laminates  at  elevated  temperatures.  They  introduced  shift  factors  for  both  the 
static  strength  and  the  fatigue  function  to  account  for  the  temperature  effect.  By 
using  lamination  theory,  fatigue  function  and  temperature  shift  factors,  the  fatigue 
lives  at  different  temperatures  were  predicted  and  correlated.  These  results  may 
permit  the  life  prediction  for  long-time  and  low  temperature  from  data  obtained  in 
short  time  at  high-temperature.  However,  all  their  work  focused  on  life 
prediction.  They  did  not  study  damage  growth,  such  as  transverse  cracking. 

As  mentioned  before,  damage  development  in  fiber  reinforced  polymer 
composite  laminates  is  a  very  complex  process.  The  static  transverse  cracking 
problem  was  studied  and  a  Monte  Carlo  simulation  and  theoretical  model  were 
developed  to  predict  crack  multiplication  with  increasing  load  at  different 
temperatures.  For  fatigue  loading,  although  the  damage  mechanisms  are  more 
complex,  experiments  show  that  the  first  stage  of  damage  still  consists  of 
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transverse  matrix  cracking.  The  formation  of  transverse  cracks  cannot  cause 
catastrophic  Mure  of  the  laminate,  but  its  presence  can  influence  the  overall 
mechanical  behavior  of  the  laminate,  such  as  stiffiiess  reduction  and  structure 
degradation.  Moreover,  the  existence  and  multiplication  of  transverse  cracks 
degrades  the  life  of  the  composite  laminate  by  inducing  other  more  severe  damage 
modes,  such  as  delamination  and  fiber  breakage. 

Experiments  show  that  crack  density  increases  with  number  of  cycles,  but 
few  works  have  been  performed  on  how  to  quantify  this  damage  growth  under 
some  specific  load  and  temperature  history.  Given  a  specific  stress  level,  stress 
ratio  and  temperature,  prediction  of  transverse  cracking  is  still  an  open  question. 
The  objective  of  the  present  work  was  to  predict  transverse  crack  development  for 
any  loading  and  temperature  history.  Since  fatigue  is  affected  by  statistical 
variability,  the  Monte  Carlo  technique  was  introduced  to  simulate  the  fatigue 
transverse  cracking  process  based  on  the  stress  analysis  model,  stress-life  curve 
and  damage  accumulation  model.  Simulation  results  were  compared  with 
experimental  ones  at  different  stress  levels  and  temperatures  and  good  agreement 
was  obtained.  Finally,  the  temperature  effect  was  highlighted  in  this  fatigue 
damage  model. 

Experimental  Procedure 

The  material  in  this  investigation  was  IM7/977-3  carbon/epoxy,  as  before. 
Unidirectional  [90s]  laminates  were  made  to  obtain  the  stress-life  curve  of  the 
unidirectional  lamina.  Two  cross-ply  panels,  [0/904]  and  [0/902]s ,  were  made  to 
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run  fatigue  tests.  All  panels  were  inspected  by  ultrasonic  C-scanning  and  found  to 
have  no  discernible  defects.  The  fiber  volume  ratio  was  0.65,  as  determined  by 
image  analysis  of  photomicrographs  of  transverse  (to  the  fibers)  cross  sections  of 
the  composite  [2].  The  [908]  and  crossply  laminate  panels  were  cut  with  a 
diamond  saw  into  22.9  cm.  (9.0  inch)  x  2.54  cm  (1.0  inch)  coupons  for  fatigue 
testing. 

To  prevent  fiber  breakage  during  gripping,  end  tabs  were  used.  For 
fatigue  testing,  especially  high  temperature  testing,  aluminum  sheets  were  used  as 
tabs  because  glass-epoxy  tabs  exhibited  cohesive  failures  during  the  fatigue  tests. 
The  choice  of  adhesive  bonding  the  tabs  to  the  specimen  is  also  very  critical  for 
high  temperature  fatigue  tests.  Experiments  showed  that  FM355  high  temperature 
adhesive  film  (Cytec  Fiberrite  Inc.)  works  well  for  high  temperature  fatigue  tests 
up  to  149°  C  (300°  F). 

All  tests  were  performed  in  an  Instron  1331  hydraulic  test  machine  with  an 
Instron  8500  controller.  An  environmental  chamber  attached  to  the  testing 
machine  was  used  for  high  temperature  tests. 

The  first  stage  of  fatigue  damage  in  a  cross-ply  laminate  consists  of 
development  and  multiplication  of  transverse  matrix  cracks.  This  process 
depends  on  the  current  state  of  stress  and  fatigue  characteristics  of  the  90°  lamina. 
The  stress-life  curve  of  the  90°  lamina  is  essential  in  characterizing  the  fatigue 
damage  behavior  of  cross-ply  laminates. 

Unidirectional  [908]  specimens  were  used  to  measure  fatigue  life  (number 
of  cycles  to  failure).  Fatigue  tests  were  conducted  using  a  sinusoidal  tension- 
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tension  cyclic  load  with  a  stress  ratio  of  R  =  0.5  at  a  frequency  of  4  Hz.  Each 
specimen  was  carefully  aligned  before  loading  to  avoid  undesirable  bending 
failure.  Both  room,  24°  C  (75°  F),  and  high  temperature,  149°  C  (300°  F),  fatigue 
tests  were  run  at  maximum  stresses  between  0.9  F2t  and  0.6  F2t.  For  each  stress 
level,  three  to  five  specimens  were  tested  to  failure  and  the  mean  value  was 
chosen  as  the  life  at  this  stress  level.  All  fatigue  tests  run  at  a  stress  of  0.6  F2t  ran 
out  at  106  cycles. 

The  fatigue  data  obtained  are  shown  in  Figure  50  and  fitted  by  a  straight 
line  for  stress  levels  above  0.7  F2t.  An  exponential  function  was  used  to  fit  the 
fatigue  data  between  0.7  F2t  and  0.6  l^t-  From  these  data  it  was  evident  that  at  the 
same  normalized  stress  level,  the  life  at  high  temperature  is  shorter  than  that  at 
room  temperature. 


Fig.  50.  Stress-Life  Curves  of  [90g]  Lamina  at  Two  Temperatures 
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After  obtaining  the  stress-life  curve  of  the  unidirectional  transverse 
lamina,  fatigue  damage  mechanisms  of  cross-ply  laminates  were  studied  both  at 
room  and  elevated  temperatures.  Two  layups  were  used,  [0/904]s  and  [0/902]s. 
One  of  the  specimen  edges  was  polished  to  12.5  microns  to  facilitate  damage 
inspection  and  crack  counting  with  a  microscope.  To  prevent  fiber  breakage 
during  gripping  and  slipping  between  specimen  and  grips,  0.81  and  0.41  mm 
(0.032  and  0.016  in.)  thick  aluminum  sheets  cut  into  3.81  x  2.54  cm  (1.5  x  1  in.) 
coupons  were  used  as  end  tabs  for  the  [0/904]  and  [0/902]s  specimens, 
respectively.  One  end  of  the  tab  was  tapered  to  reduce  stress  concentration.  Tabs 
were  bonded  to  the  specimens  with  FM355  adhesive  film  (Cytec  Fiberite  Inc.), 
which  is  highly  recommended  for  fatigue  testing.  To  prevent  cohesive  failure 
between  tabs  and  adhesive,  the  surface  of  the  end  tabs  in  contact  with  the 
adhesive  was  cleaned  thoroughly  with  MEK.  In  addition,  to  enhance  bonding 
between  specimen  and  adhesive  film,  an  approximately  5  cm  (2  in.)  long  section 
at  each  end  of  the  specimen  was  roughened  with  300  grit  sand  paper.  The 
adhesive  film  was  cured  at  177°  C  (350°  F)  for  one  hour  before  testing. 

All  fatigue  tests  were  performed  in  an  Instron  1331  hydraulic  testing 
machine  with  an  Instron  8500  controller  using  the  same  stress  ratio,  R=0.5,  in  the 
transverse  layer  as  in  the  case  of  the  unidirectional  lamina  tests,  at  a  frequency  of 
4  Hz.  Three  specimens  were  tested  under  each  test  condition.  Tension-tension 
fatigue  tests  were  performed  at  two  local  stress  levels  for  each  layup  at  different 
temperatures:  G2xmaX  /F2t  =  1  *1  and  o,2xmax  /  F2t  =  1  -4  lbr  [0/904]s  and 
a2xmax/F2t  =1.2  for  [0/902]s,  where  a2xmax  is  the  maximum  applied  stress  in 
the  transverse  layer  and  F2t  is  the  static  transverse  tensile  strength.  High 
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temperature  tests  were  conducted  in  an  environmental  chamber  attached  to  the 
testing  machine. 

Fatigue  tests  were  interrupted  after  the  first  cycle  to  determine  initial  crack 
density  and  after  some  predetermined  number  cycles  to  monitor  crack 
development  and  multiplication.  Then,  the  specimen  was  taken  to  an  optical 
microscope  to  observe  the  damage  from  the  edge.  Cracks  in  the  90°  layer  were 
counted  over  the  gage  length.  It  was  found  that  almost  all  cracks  extended 
through  the  thickness  of  the  90°  layer,  and  new  through  the  thickness  cracks 
appeared  instantaneously  without  a  noticeable  crack  propagation  stage.  Figure  51 
shows  typical  images  of  cracking  in  the  90°  layer  of  [0/904]s  laminates  by  X- 
radiography  at  various  stress  levels  and  temperatures.  X-radiography  shows  that 
all  fatigue  cracks  extended  through  the  width  of  the  specimen.  Counting  cracks 
on  the  edge  was  much  simpler  than  using  other  techniques,  such  as  replication  and 
X-radiography  to  obtain  crack  density.  Each  fatigue  test  was  run  for  up  to  106 
cycles  and  crack  density  at  every  predetermined  number  of  cycles  was  recorded. 
During  the  test,  no  delamination  at  interfaces  was  noticed. 

Damage  Development 

Characteristics  of  Fatigue  Cracking  in  Crossply  Laminates 

The  measured  crack  density  in  the  transverse  layer  was  plotted  versus 
number  of  cycles  in  Figs.  52-57,  for  the  two  laminates  tested,  at  three  stress  levels 
and  two  temperatures.  Each  point  on  these  curves  represents  an  average  from 
three  tests. 
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Fig.  51.  X-Radiographs  of  [0/904]s  IM7/977-3  Carbon/Epoxy  Laminate  under 
Fatigue  Loading  at  Various  Applied  Stress  Levels  and  Temperatures 
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Fig.  52.  Crack  Density  vs.  Number  of  Cycles  for  [0/904]s  Laminate  at 
a™*  /  F2t  =  1 .1  and  24°  C  (75°  F) 


Fig.  53.  Crack  Density  vs.  Number  of  Cycles  for  [0/904]s  Laminate  at 
arx/F2t  =  1.4  and  24°  C  (75°  F) 
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54.  Crack  Density  vs.  Number  of  Cycles  for  [0/904]s  Laminate  at 
/p2t  =  U  and  149°  C  (300°  F) 


55.  Crack  Density  vs.  Number  of  Cycles  for  [0/904]s  Laminate  at 
a”3* /Fit  =  1.4  and  149° C (300° F) 


56.  Crack  Density  vs.  Number  of  Cycles  for  [0/902]s  Laminate  at 
og*  /  F2t  =  1 .2  and  24°  C  (75°  F) 


57.  Crack  Density  vs.  Number  of  Cycles  for  [0/902]s  Laminate  at 
<52T  /F2t  =  1.4  and  24°  C  (75°  F) 
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From  the  experimental  results,  it  was  observed  that  the  crack  density  rate 
with  number  of  cycles  is  not  constant  during  the  cracking  process.  Generally, 
there  were  three  damage  stages.  The  first  one  corresponds  to  cracking  of 
segments  whose  length  (crack  spacing)  is  well  above  the  average  crack  spacing 
seen  after  the  first  cycle.  At  this  stage,  the  increase  in  crack  density  (crack 
number)  with  number  of  cycles  is  slow.  The  characteristic  of  the  second  damage 
stage  is  a  steep  increase  in  crack  density,  which  is  attributed  to  multiple  fatigue 
cracking  of  segments  of  length  approximately  equal  to  the  average  crack  spacing 
after  the  first  cycle.  In  the  third  stage,  the  crack  density  approaches  a  saturation 
limit  and  the  slope  of  the  curve  begins  to  level  off.  In  fatigue  loading,  it  is 
believed  that  crack  saturation  is  related  to  the  endurance  limit  of  the  stress-life 
curve  of  the  unidirectional  lamina.  Once  the  maximum  transverse  stress  in  the 
segment  falls  below  the  endurance  limit,  there  is  no  further  cracking  in  this 
segment.  It  was  also  found  out  that  the  saturation  crack  density  varies  with 
applied  stress  level. 

Another  interesting  phenomenon  studied  was  the  temperature  effect  on 
fatigue  damage.  Under  the  same  normalized  stress  level,  experimental  results 
showed  that  crack  density  at  high  temperatures  was  higher  than  that  at  room 
temperature  for  the  same  number  of  cycles.  Thus,  high  temperature  was  found  to 
be  an  accelerating  factor  in  transverse  cracking.  Since  the  damage  mechanisms 
are  the  same,  it  is  possible  to  predict  long-term,  low  temperature  fatigue  damage 
from  short  time,  high  temperature  tests.  Basically,  this  is  the  long-term  damage 
prediction  problem  which  is  very  important  for  reliability  and  durability  studies. 
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In  the  case  of  a  crossply  laminate  under  fatigue  loading,  it  is  assumed  that 
matrix  cracking  continues  in  the  transverse  layer  until  the  maximum  stress  in  the 
cracked  lamina  reaches  the  value  of  the  endurance  limit  of  the  transversely  loaded 
lamina.  The  corresponding  crack  density  at  this  characteristic  damage  state 
(CDS),  Xcds,  can  be  determined  from  the  stress/crack  density  relation  (Eq.  24), 
where  the  lamina  stress  is  replaced  by  the  endurance  limit  stress,  ce/. 
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The  above  relation  can  be  normalized  by  the  corresponding  quasi-static  transverse 
tensile  strength  F2t.  It  is  also  noted  that  in  the  above  equation  the  average 
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The  normalized  endurance  limit  is  obtained  or  estimated  from  the  S-N 
curve  of  Fig.  50  and  used  with  eq.  (31)  to  determine  the  corresponding  crack 
density, 
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The  corresponding  number  of  cycles,  Ncds?  can  be  obtained  or  estimated 
from  the  S-N  curve  of  the  unidirectional  lamina  (Fig.  50). 


72 


Experimental  measurements  of  crack  density  and  number  of  cycles  are 


plotted  in  normalized  form  as 

X-Xn  logN 

^cds  “  logNCDS 

and  shown  in  Fig.  58,  where  X0  is  the  initial  (first  cycle)  crack  density,  Xcds  is  the 
calculated  value  from  eq.  (32).  In  the  present  case  the  normalized  endurance  limit 
used  was 


^sL  =  0.60 

ht 

Results  for  different  temperatures  and  stress  amplitudes  appear  to  be  reasonably 
represented  by  a  master  curve  as  shown  in  Fig.  58.  Furthermore,  it  appears  that 
results  for  elevated  temperatures  correspond  to  a  later  stage  in  fatigue  life,  a  result 
of  potential  use  for  accelerated  testing  and  life  prediction  methodologies. 


Fig.  58  Predicted  and  Measured  Normalized  Transverse  Crack  Density  versus 
Normalized  Life  in  Crossply  Laminate  under  Tension-Tension  Fatigue 
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Monte  Carlo  Simulation  Approach 

In  order  to  develop  and  apply  a  statistical  model  to  the  onset  and 
multiplication  of  transverse  cracks  the  following  are  needed:  (a)  criteria  for 
transverse  cracking  under  fatigue,  (b)  a  damage  accumulation  theory  or  rule  and 
(c)  a  stress  analysis  model. 

Currently,  there  are  two  approaches  explaining  fatigue  cracking.  In  one 
approach,  based  on  fracture  mechanics,  it  is  assumed  that  new  cracks  start  from 
some  effective  pre-existing  flaws  which  increase  in  size  with  number  of  cycles 
[27].  When  a  flaw  reaches  a  characteristic  size,  a  new  transverse  crack  is  formed. 
Flaw  growth  during  fatigue  is  described  by  Paris’  law.  However,  this  approach 
requires  an  accurate  determination  of  fracture  toughness  parameters  which  are  not 
easy  to  obtain  and  thus,  resort  is  made  to  empirical  factors.  Furthermore,  this 
approach  focuses  mainly  on  the  crack  propagation  stage.  Experimental 
observations  show  that  transverse  cracks  propagate  across  the  thickness  and  width 
nearly  instantaneously,  thus,  minimizing  the  importance  of  the  propagation  stage. 
The  alternative  approach  is  based  on  the  S-N  curve  which  relates  failure  in  the  90° 
layer  (cracking)  to  the  state  of  stress  and  number  of  cycles. 

In  order  to  account  for  the  changing  state  of  stress  (due  to  progressive 
cracking)  and  the  varying  damage  history  at  every  point  of  the  90°  layer  in  the 
specimen,  there  is  need  for  a  stress  analysis  model  and  a  damage  accumulation 
model. 

The  statistical  model  proposed  in  this  study  is  based  on  the  S-N  curve  of 
the  unidirectional  90°  lamina  (cracking  criterion),  a  damage  accumulation  model 
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(Miner’s  rule)  and  a  stress  analysis  model  based  on  the  shear  lag  analysis 
discussed  before  [22].  Monte  Carlo  simulations  were  performed  for  different 
crossply  laminates,  load  levels  and  temperatures. 

The  stress  analysis  model  discussed  before  (eqns.  16  and  17)  was  used  to 
determine  the  cyclic  stresses  in  the  cracked  90°  layer.  The  maximum  cyclic  stress 
in  the  90°  layer  between  two  cracks  is  given  by 
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where  ct  and  a  are  the  maximum  applied  cyclic  stresses  in  the  laminate 

and  90°  layer,  respectively.  If  a^1xax>  F2t  ,  there  will  be  some  initial  cracking 

after  the  first  cycle.  The  following  discussion  focuses  on  this  case. 

A  basic  problem  in  fatigue  damage  is  quantification  of  cumulative  damage 
and  life  prediction  under  variable  amplitude  cyclic  loading.  For  a  composite 
laminate,  such  as  a  crossply  laminate,  even  if  the  global  applied  fatigue  load  level 
is  fixed,  the  load  share  of  each  layer  changes  continuously  due  to  crack 
multiplication  in  the  transverse  layer,  which  also  changes  the  local  axial  stresses 
in  the  transverse  layer.  Thus,  the  effect  of  varying  stress  level  on  the  fatigue 
damage  accumulation  has  to  be  addressed  in  order  to  predict  the  transverse 
cracking  progress.  A  simple  and  well-known  damage  accumulation  model  is 
Miner’s  rule  [47].  According  to  this  rule,  the  cumulative  damage  at  a  point  in  the 
specimen  subjected  to  k  blocks  of  fatigue  loading  at  various  stress  levels,  Oi, 
a2, . ,  ork,  for  corresponding  cycles,  ni,  n2, . %  is 
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where  Ni  is  the  number  of  cycles  to  failure  at  the  point  under  stress  Oj.  A  crack 
will  be  produced  at  the  point  of  interest  when 
Dk  =  l 

Miner’s  rule  does  not  account  for  the  effect  of  loading  sequence.  Other 
nonlinear  fatigue  models  considering  load  sequence  effects  have  been  proposed 
[33,  36,  38,  44,  45].  Hashin  reviewed  the  cumulative  damage  problem  and 
compared  predictions  of  nonlinear  models,  including  his  own,  with  experimental 
results  [48],  He  found  that  none  of  the  above  models  gave  satisfactory  predictions 
in  all  cases.  In  particular,  none  of  the  nonlinear  models  gave  a  better  fit  than 
Miner’s  rule  in  every  case.  Actually,  all  the  prediction  differences  are 
overshadowed  by  the  fatigue  data  scatter.  Thus,  Miner’s  rule  was  selected  as  an 
adequate  fatigue  damage  accumulation  model  in  the  present  study. 

The  first  set  of  cracks,  not  necessarily  equidistant,  at  the  first  fatigue  cycle 
was  determined  by  the  simulation  method  used  for  the  static  loading  case.  The 
initial  crack  density,  X0,  and  the  location  of  each  crack  were  obtained. 

The  first  loading  block  at  an  applied  laminate  stress  ax  covers  the  period 
up  to  the  formation  of  the  next  crack.  The  critical  point  in  each  segment  of  the 
90°  layer  between  cracks  is  identified  and  the  numberof  cycles  needed  to  produce 
a  crack  at  that  point  is  determined  from  the  S-N  curve  for  the  lamina  (Fig.  50). 
The  next  crack  will  appear  at  the  point  needing  the  fewest  number  of  cycles  to 


failure. 


In  the  second  loading  block,  up  to  the  formation  of  the  next  crack,  the 
cumulative  damage  at  the  critical  points  of  the  uncracked  segments  will  continue 
to  be  monitored  to  determine  the  number  of  cycles  needed  to  form  a  crack. 
Simultaneously,  all  points  in  the  newly  formed  segments  (on  either  side  of  the 
new  crack)  will  be  interrogated  to  determine  the  critical  points  within  them 
(minimum  remaining  life).  The  next  crack  will  appear  at  the  critical  point  with 
the  shortest  remaining  life  (smallest  number  of  cycles)  among  the  set  of  critical 
points. 

In  general,  at  a  given  loading  block  k,  the  number  and  location  of  cracks 
are  known.  At  the  end  of  the  previous  loading  stage,  or  block  k-1,  the  crack 

density  Xk_,  and  the  accumulated  damage  D^_,  are  known  at  every  pointy.  The 
remaining  life  (number  of  cycles)  during  loading  stage  k  at  every  critical  pointy  is 
obtained  by  calculating  the  local  stress  (<*2x)i  from  the  shear  lag  model,  S-N 

curve  and  damage  accumulation  model.  The  remaining  life  n[  at  the  critical 
point  j  is 

^-=l-Dti_,  (35) 

NJ 

where  N  (  and  N  are  lifetimes  obtained  from  the  S-N  curve. 

The  new  crack,  corresponding  to  crack  density  X*,  will  form  at  the  point  r 
with  the  minimum  n  k  remaining  life  during  loading  stage  k.  The  cumulative 
number  of  cycles  corresponding  to  A*  is 

NcL  =  Nc^1+min(nJk)=NcL1+nL  (36) 
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and  the  accumulated  damage  at  a  critical  point  j  in  stage  k  is 

D{=^  +  D,j_,  (37) 

Nk 

This  process  can  be  continued  up  to  the  point  when 
^  =  ^cds 

and  Ncum  =  NCDS 

Transverse  cracking  of  two  cross-ply  laminates,  [0/904]s  and  [0/902]s,  for 
two  different  load  levels  and  two  different  temperatures,  was  simulated. 
Simulation  results  were  compared  with  experimental  results  as  shown  in  Figs.  51- 
56.  The  simulation  results  using  Miner’s  rule  agree  well  with  the  experimental 
data. 

SUMMARY  AND  CONCLUSIONS 

An  investigation  was  conducted  on  characterization  of  constitutive 
behavior,  physical  aging,  failure  mechanisms  and  damage  accumulation  in  high 
temperature  polymer  matrix  composites. 

Most  aging  studies  deal  with  aging  effects  on  material  properties,  such  as 
creep  and  usually  under  idealized  isothermal  conditions.  It  is  difficult  to  predict 
long  term  behavior  of  composites  undergoing  loading  and  environmental 
fluctuations.  Instead  of  the  effects,  the  aging-induced  physical  changes 
(volumetric  shrinkage),  which  describe  the  “degree”  of  aging  and  not  just 
chronological  age,  were  measured  in  this  study.  A  physical  aging  model  was 
proposed  that  can  describe/predict  aging  under  any  temperature  history. 
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Furthermore,  a  new  dynamic  pulse  loading  method  was  developed  for  assessing 
aging-induced  viscoelastic  effects  and  relating  them  to  volumetric  shrinkage. 

A  great  deal  of  work  has  been  reported  on  modeling  failure  and  damage 
development  in  composites.  The  eternal  challenge  is  to  predict  long  term 
behavior  based  on  short  term  evidence.  A  unique  and  simple  general  relationship 
was  developed  between  applied  load  and  crack  density  in  composite  laminates, 
independent  of  lamination  geometry  and  temperature  and  accounting  for 
statistical  variability.  Under  fatigue  conditions  a  master  curve  was  developed 
relating  damage  (crack  density)  and  fatigue  cycles,  independent  of  layup  and 
temperature.  These  results  will  help  in  the  development  of  accelerated  test 
methods  and  life  prediction  models. 

The  research  results  obtained  to  date  can  serve  to  develop  design 
methodology  for  polymer  composite  structures  for  long  term  service  under 
varying  environmental  conditions.  In  cases  where  physical  aging  is  a  significant 
factor,  there  is  a  way  to  assess  the  current  degree  of  aging,  irrespective  of  and 
without  knowing  the  past  history.  Accelerated  test  methodology  could  be  based 
on  understanding  the  effects  of  temperature  and  load  amplitude  on  fatigue  life  and 
failure  mechanisms. 
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